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Abstract 

We show that atoms from wide velocity interval can be concurrently involved in 
Doppler-free two-photon resonant far from frequency degenerate four-wave mixing with 
the aid of auxiliary electromagnetic field. This gives rise to substantial enhancement of the 
output radiation generated in optically thick medium. Numerical illustrations addressed 
to typical experimental conditions are given. 

PACS number(s): 42.50. Gy, 42.65.Dr, 42.65.Ky 

Strong optical resonances inherent to free atoms and molecules are negated by the fact that 
only small fraction of the species can be concurrently resonance coupled in warm bulk gases. 
This is because of Maxwell distribution of the Doppler shifts of their resonances. Doppler-free 
(DF) coupling is usually achievable under equal frequency counter-propagating weak waves 
in two-photon-resonant ladder schemes. Because of the phase matching requirements, such 
schematic can not be implemented for far from frequency degenerated four- wave mixing (FWM). 
Second, the conditions of intermediate one-photon quasi-resonance can be very seldom satisfied 
in this case. Third, DF coupling vanishes with the growth of not only frequency difference but 
with strengths of the coupled fields too due to ac-Stark effect. Forth, DF coupling can not be 
routinely achieved in Raman schemes, like considered in this paper, because of the inherent 
frequency difference. Fifth, DF absorption does not indicate readily achievable DF FWM 
polarization in general case. Moreover, it is not obvious that increased FWM polarization 
would result in enhanced output of generated radiation, because of accompanying increased 
absorption. 

This paper is aimed at demonstrating that the outlined limitations in optical physics can 
be removed by making use quantum coherence processes induced by an auxiliary intense con- 
trol field. Considered effects lead to concurrent contribution of atoms from a wide velocity 
interval to the induced resonance and to eliminating it's Doppler broadening under moderate 
light intensities. We investigate two- photon-resonant Raman-type FWM, controlled by an 




Figure 1: Coupling scheme ((a)) and enhancement of quantum conversion efficiency in opti- 
cally thick medium with the aid of counter-propagating control field ((b) and (c) ). (ao\Z - 
optical thickness of the medium at frequency u;oi under E 2 = E 3 (oj 3 ) = 0). (b): G\ 2 = 128.5 
MHz, G 23 = 25.2 GHz, fi 2 = 92.3 GHz, fi+ = -7.3 GHz, tt 3 = 73.2 GHz; inset: same but 
E- = 0. (c): G l2 = 74.2 MHz, G 23 = 635.8 MHz, fi 2 = 2.3 GHz, fi+ = -1.96 GHz, fij = 1.83 
GHz; plot b: same but E 3 = 0; inset: G 12 = 742.2 MHz, G% 3 = 0, fi 2 = 2.3 GHz, fi+ = 0. 
In all cases G^ 3 = 5.78 MHz, oj\ is set to induced Doppler-free resonance (for insets - to the 
maximum of nonlinear susceptibility). 

auxiliary driving field. The explicit formulae for power- dependent absorption/gain indices and 
for nonlinear FWM susceptibilities, accounting for interplay of power and Doppler shifts of 
the resonances and illustrating the major idea of the proposed method are derived. In order 
to avoid accompanying population transfer, which would complicate the formulae and would 
mask the major effect under consideration, the field coupled to the ground state is assumed to 
be weak. 

A possible achievement of sub-Doppler resolution using intense control field was shown in 
U, [2|, [3], f|, |5j. We propose a novel scheme, which enables to control FWM coupling with 
the aid of auxiliary electromagnetic (EM) field, taking no part in the FWM process itself. 
Accompanying increase of absorption of the fundamental radiation is considered too. As the 
outcome, substantial enhancement in quantum conversion efficiency in optically thick Doppler- 
broadened medium is shown. Numerical illustrations are given for the model, relevant to the 
FWM experiments with sodium dimer vapors 0. As a matter of fact that detuning from the 
intermediate resonance are larger than the Doppler width of the transition and the populated 
ground level is coupled to weak fields only, none of the CPT or EIT effects, usually employed 
for the enhancement of resonant FWM J/J , are involved in the proposed technique. 

Consider FWM process uj\ — uj 2 + UJ 3 =^5 = ^4 and transition configuration (fig. la), 
similar to those studied in the experiments ||. However, in our case the EM radiation E 3 (t, z) 
consists of two components: weak £"^(0;^") and counter propagating strong one E 3 {uj 3 ). Their 
frequencies can be the same or different. Strong radiation E 2 and weak E\ co-propagate in the 
same direction as E 3 . Only lower level remains populated, because Ei is assumed so weak, that 
the populations can not be driven. Density matrix equations in the interaction representation 



are: 



AnPoi = i{pooVbi + P02V21}, (1) 
^03P03 = i{pooV(Q + P02(V£ + V 23 )}, (2) 
L 2p02 = i{p0lVl2 + PtaOQ + V32)}, (3) 

where = d/dt + vV + T^, = exp{z(f2jt — fcjz)}, V^g = exp{i(Q,ft =F fcj^)}) 
= —Ejdij/2h, G 23 = —E^d 23 /2h - are coupling Rabi frequencies, fij - are corresponding 
resonance detunings (e.g., fii = u± — Woi), - homogeneous half- widths of the transitions. 

As follows from (|l|) - @, induced atomic coherence P02 gives rise to the components in 
polarizations, responsible for novel effects in absorption and generation of the radiations under 
consideration. In the lowest order on the strength of the weak fields solution of the equations 
(HP ~~ (0) can be found in the form: 

P02 = ?02 exp{i [(Qi - Q 2 )t — (fci — k 2 )z] } + r^ 2 exp{i[(Q 4 — fig )t — (k 4 — k 3 )z] } + Tq 2 exp{i - 
^3 )t~ (ks + K) z ]}i P03 = r 03 exp{i(tt 4 t - k^z) + f 03 exp{i(D, 4 t- k s z)} + fQ 3 exp{i(^t- k^z)}, 
P01 — r oi exp{i(Qi — kiz)}+f i exp{i(r2xt — A;i2)}+f 01 expj^fift — fcf z)}, where f2j = fli — Q 2 + 
fig , = ^4 — ^3 +^2, &s = ki — k 2 + k 3 , kg = k\ — k 2 — k 3 ,k\ = k 4 — k 3 +k 2 ,k± = k 4 + k 3 +k 2 . 

Equations for the density-matrix amplitudes become algebraic. With aid of solution for 
^035 r oi and ?"o3 expressions for the susceptibilities, dressed by the strong fields E 2 and E 3 , can 
be routinely obtained and presented as: 
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where X01 and X04 are corresponding resonant values under the strong fields being turned off, 

P01 = r i + - kxv),p^ = r i + i(^r - Kv),p 03 = r 03 + ayt 4 - M),^! = r o3 + 

i(fl 4 - k 4 v), P 02 = r 02 + - O2 - - fc 2 )u], ^02 = r 02 + i[^4 -^3 ~ (^4 + fcg )u], P 2 = 
-P02 + \Gi 2 \ 2 /Poi + \G 23 \ 2 / P^, v is projection of atom velocity on z. Difference between k\ and 
k\ as well as between k 4 and fcs is neglected here. 

With account of absorption but neglecting depletion of fundamental radiations due to FWM 
conversion, reduced equation for E 4 can be written as: 

dE 4 (z)/dz = i2nhxf ) E 1 (0)E;E 3 exp(-iAKz), (7) 

where AK = K 4 — K\ + K 2 — K 3 , Kj = kj — iaj/2 - are complex wave numbers, ctj - power- 
dependent absorption indices. Quantum conversion efficiency (QCE) of E\ into E 4 along the 
medium rj q (z) is given by the expression: r] q (z) = (ui/u 4 )\E 4 (z)/ Ei(0)\ 2 exp(— a 4 z). 
From (0) one obtains: 

= itr^ exp(-a 4 z)|exp(-iAAT2:) - l| 2 . (8) 

CU4 Mr ' ' 



Pre-exponential factor can be expressed over Rabi frequencies, reduced nonlinear susceptibility 
and absorption indices considered below, ratios of the transition widths and | g?03 1 2 / 1 <^oi 1 2 - The 
last factor is proportional to the ratio of the spontaneous relaxation rates. Thus QCE can be 
found as absolute value, dependent on the optical thickness of the medium. 

The major physics underlying the proposed technique is as follows. Modulation of the atomic 
wave-functions by the driving fields gives rise to the Autler-Townes splitting, which exhibits 
itself in our case as resonance shift. Besides intensities, the later depends on detunings of the 
driving fields and consequently - on the atomic velocities. It turns out that under appropriate 
intensities the resonances of atoms at different velocities can be shifted to the approximately 
same position. To illustrate that, consider one-photon detunings, substantially greater than 
corresponding Doppler HWHM. Then the resonance power-shift factors in (|j) - can be 
presented as: \G\ 2 /P ~ (1 + ikv/p)\G\ 2 /p, where p is corresponding factor P at v — 0. This 
shows possible control of the resonance Doppler shifts through the power shifts. More details 
can be found in 0, m the same way a factor in the denominators of (§]), @, indicating 
dressed two-photon resonance, can be presented as: 

P 02 ~ f 02 + ifl 02 - t{(\G 12 \ 2 /Q 2 1 )k 1 + 

where r 2 and ^02 give half- width and position of the induced resonance. As follows from (0), 
under proper choice of detuning, relative propagation direction and intensity of the control field 
E% , all Doppler shifts can be compensated by the power shifts in a such way, that dependence 
on v vanishes in the given linear on v approximation. This indicates trapping of all atoms, 
independent of their velocities in DF dressed two-photon resonance. It is seen, that as a 
matter of fact that k 2 < k±, elimination of Doppler broadening is not possible in the schematic 
under consideration with the aid of only driving field E 2 . However it becomes possible with 
an auxiliary counter-propagating control field E%, which does not contribute directly in FWM 
because of phase mismatch. The equations (f|), (|^) show similar behavior of absorption index 
and nonlinear susceptibility near induced resonance. 

While approaching closer to the intermediate resonances, required intensities become lower, 
but relative contribution of the neglected terms, proportional to the higher orders on kiv/Qi, 
grows. This leads to decrease of the coherently coupled velocity interval. 

The discussed outcomes can be illustrated with the numerical model of sodium dimer tran- 
sitions 0: A01 = 661 nm, A 12 = 746 nm, A23 = 514 nm and A03 = 473 nm. Corresponding 
homogeneous half-widths of the transitions are 20.69, 23.08, 18.30 and 15.92 MHz, Doppler 
HWHM - 0.678, 0.601, 0.873 and 0.948 GHz. 

Figure 2 depicts contribution of molecules at different velocities to the absorption index 
oti(uJi) ~ -Re{xi/xoi} an d to nonlinear susceptibility (trivial Maxwell envelopes are removed), 
while conditions of elimination of Doppler broadening are fulfilled. The figure shows potentials 
of coherent coupling of molecules from wide velocity interval compared with the width of the 




Figure 2: Velocity distribution of the squared modulus of FWM nonlinear susceptibility ((a)) 
and absorption index at U\ ((c)) in Doppler-free resonance; (b) and (d) - corresponding 
velocity-averaged sub-Doppler resonances (scaled to the corresponding value at the frequency 
of control-field induced resonance but under E^(u^) =0). u - thermal velocity, Aum - 
Doppler HWHM of the transition 01. The insets - same functions at E^(lj^) = 0. Detunings 
and intensities are the same as in fig. lb. 

Maxwell distribution, unlike the case in the absence of the control field. This gives rise to 
strong sub-Doppler resonances. Figure 3 shows modification of the effects while tuning closer 
to the intermediate resonances. Despite the growth of absorption a\, the proposed manipulating 
results in substantial increase in output of generated radiation at us (fig-1 (b) and (c)). 

In conclusion, we show that substantial enhancement in nonlinear-optical response of a 
Doppler broadened medium can be achieved by coherent driving of quantum transitions so 
that molecules from wide velocity interval become trapped to one and the same dressed two- 
photon Doppler-free resonance. The required intensities can be decreased by tuning driving 
frequencies closer to one-photon resonances, while the coupled velocity interval decreases too. 
For ladder-type schemes, where Doppler-broadening of two-photon resonances is much larger 
compared to Raman-like schemes, the considered effects are even more pronounced. 
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Figure 3: Velocity distribution of the squared modulus of FWM nonlinear susceptibility ((a)) 
and absorption index at ui ((c)) in Doppler-free resonance; (b) and (d) - corresponding 
velocity-averaged sub-Doppler resonances (scaled to the corresponding value at the frequency of 
control-field induced resonance but under E^(u^) = 0). The insets: same profiles at E^(u^) = 
0. Detunings and intensities are the same as in fig. lc. 
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